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We consider variational (density functional) models of fluids confined in parallel-plate geometries 
(with walls situated in the planes z = and z — L respectively) and focus on the structure of the 
pair correlation function G(ri,r2). We show that for local variational models there exist two non- 
trivial identities relating both the transverse Fourier transform G(z M , z„; q) and the zeroth moment 
Go(Z(„ z v ) at different positions Zi, 22 and 23. These relations form an algebra which severely restricts 
the possible form of the function Go(z l _ l , z v ). For the common situations in which the equilibrium 
one-body (magnetization/number density) profile mo (z) exhibits an odd or even reflection symmetry 
in the z = h plane the algebra simplifies considerably and is used to relate the correlation function 
to the finite-size excess free-energy j(L). We rederive non-trivial scaling expressions for the finite- 
size contribution to the free-energy at bulk criticality and for systems where large scale interfacial 
fluctuations are present. Extensions to non-planar geometries are also considered. 



I. INTRODUCTION 

In this paper we discuss the structure of correlation 
functions for fluids (or Ising-like magnets) adsorbed at 
walls and confined in parallel-plate (thin film) , cylindrical 
and spherical geometries. Such systems have been 
extensively studied in recent years to assess the influence 
that surface and finite-size effects have on bulk phase 
coexistence and criticality. Here we are not primarily 
interested in the phenomenology but seek rather to 
understand whether a central assumption crucial to a 
broad class of widely used variational (density functional) 
models has any consequences for the behaviour of the 
correlation function G(ri,r 2 ). Specifically, we shall 
show that for local functional models there exist very 
general relations which may be viewed as an algebra 
restricting the possible structure of G(ri,r 2 ). These 
turn out to be particularly powerful when the one- 
body (magnetization/density) profile has an odd/even 
reflection symmetry (as is often the case) and enable 
us to re-derive non-trivial scaling laws for the finite-size 
contribution to the free-energy for systems close to the 
critical point [fi],|| or that have large scale interfacial 
fluctuations H|§. 

The starting point of our analysis is an appropriate 
variational model for the free-energy. For almost all 
of our article we will consider grand potential-like 
functionals f2[m(r)], where m(r) denotes the appropriate 
order parameter or local density variable. Thus for 
fluid systems m(r) corresponds to the number density 
although we shall adopt a magnetic notation and will 
often refer to m(r) as the local magnetization or spin 
density. The grand potential functional is written 



n[m(r)] = F[i 



i{v)]- J>, 



i(r)h(r)dr 



(1.1) 



where J-[m(r)] is the Helmholtz free-energy functional 
and h(r) is the local magnetic field at point r. In a fluid 



context h(r) may be identified with /1 — V^ x t(r) where 
fi is the chemical potential and V^ x t(r) is the external 
field generated by the confining walls. Minimization 
of Q[m(r)] with respect to magnetization configurations 
recovers the thermodynamic grand potential fl ^-0] . 



O = minf2[m(r)] 
= n[m Q (r)] 



(1.2) 



where mo(r) is the equilibrium magnetization profile. 

The connected correlation function G(ri, r 2 ) is defined 
in the usual way by 

G(ri,r 2 ) = (m(ri)m(ri)) - (m(ri ))(m(r 2 )) (1.3) 

and satisfies the Ornstein-Zernike integral equation jB) 



J dv'C{v 1 ,v')G{v',v 2 ) = 6{v 1 -v 2 ) 



(1.4) 



where the direct correlation function G(ri,r 2 ) is defined 
by 



G(r 1; r 2 ) = 



1 



S 2 F[m{r)} 



ksT 5m(ri)5m(r 2 ) 



(1.5) 



and is evaluated at equilibrium. 

One of the chief merits of the variational approach is 
the prescription for calculating G(ri,r2) via the direct 
correlation function route embodied in (1.4) and ( |l.5| ). 
Often the fluid has a symmetry which makes the task 
simpler. For example, in planar parallel-plate geometries 
the one-body profile is a function of one coordinate, 2 
say, and it is convenient to Fourier transform G(i"i,r 2 ) 
and G(i"i, r 2 ) with respect to the transverse displacement 
vector yi2 of the two points ri and r 2 . Thus we define 

G(z 1 ,z 2 ;q) = J dy 12 exp(iq.y i2 )G(ri, r 2 ) (1.6) 



and C(zi,z 2 ;q) similarly, which from (1.4) satisfy 



1 



/ dz 3 C(z 1 ,z 3 ;q)G(z 3 ,z 2 ;q) = 5(zx - z 2 ) (1.7) 

We shall be particularly interested in the zeroth moment 
corresponding to q = and define (in standard notation) 



Go(zi,z 2 ) = G(zi,z 2 ;0) 



(1.8) 



For the planar parallel-plate geometry the finite-size 
contribution to the grand potential f2 is conveniently 
measured by the surface excess quantity 



A 



(1.9) 



where V = AL is the volume of the thin film of width L 
and area A. Here u>b is the bulk grand potential density. 
The quantity 7(L) may be regarded as the finite-size 
dependent surface tension (excess free-energy per unit 
area) of the system. Differentiating with respect to L 
yields the solvation force H 



fs(L) 



dl{L) 
dL 



(1.10) 



which is a useful measure of phase behaviour in the 
confined fluid. The most important conclusions of this 
paper concern the derivation of elegant relations between 
the correlation function Gq(z\,z 2 ) and (derivatives of) 
the free-energy 7(L). 

Of course the exact density functional for non-trivial 
(interacting) three dimensional systems is not known and 
approximate models must be considered for the purposes 
of calculation. The simplest approach (and one which 
has given invaluable insights into a number of problems 
H ) is to assume that the Helmholtz functional is local so 
that we may write 



Q\m 



(r)]= j dy j dz{£ (6) (m(r),Vm(r)) 

+8(z)4>x(m(r)) + 8{z - L)(p 2 {m(r))X ( 



1.11) 



where (m, Vm) is an appropriate bulk free-energy 
density while 4>i{m) models (short ranged) interactions 
with the walls situated in the planes z — and z = L, 
say. Most often (m, Vm) is further approximated by 
a Landau expansion |6| 

C [b) (m, Vm) = -(Vm) 2 + —m 2 + —m^ - hm (1.12) 
2 2! 4! 

where, for fluid systems m must now be interpreted as 
the number density relative to the bulk critical value. 
The inadequacies of the Landau-type approach near the 
bulk critical temperature are well known though more 
general expressions for (m, Vm) can account for non- 
classical critical exponents in a phenomenological manner 
Such, generalised, local free-energy functionals have 
provided a profitable way of deriving critical exponent 



and scaling relations, although recently, it has been 
argued that local entropy-like functionals may well be 
better candidates for describing critical effects p0[ . As 
well as being restricted to systems with short ranged 
forces the local approximation is also inadequate for 
describing the pronounced oscillations in the density 
profile that occur when a high density fluid is confined in 
a narrow geometry Q . Nevertheless due to the continued 
wide spread use of local functionals we believe that 
it is worth considering in detail what restrictions the 
local approximation places on the structure of correlation 
functions. 

To address this question we borrow recent results and 
methods developed for analysing correlation functions at 
wetting transitions which have led to the introduction 
of coupled effective Hamiltonians ]Tl[— |l3|] . While we 
shall not be particularly concerned with the wetting 
transition in this article, our analysis involves the further 
development and generalization of the stiffness-matrix 
formalism present in [12| and The essence of 

the method is to consider the properties of constrained 
functionals constructed by partial minimization of the 
grand potential functional f2[m]. The Ornstein-Zernike 
equation is then recast in a convenient matrix rep- 
resentation which allows us to exploit the separable 
properties of the constrained functionals. We emphasise 
here that the partial minimisation procedure is exact in 
the present variational formalism and contrasts with the 
analogous result in the effective Hamiltonian theory of 
wetting which arises due to a saddle point approximation 
]l4| |. The present analysis also differs from the coupled 
effective Hamiltonian theory since it is necessary to 
consider the properties of constrained functionals of three 
(or more) collective coordinates rather than just two. 
For the greater part of our article, embodying Sec. |J 
~\ , we deal with planar inhomoge neous fluids described 
by local functionals of the form ( 1.11 ) but leave 
and <t>i arbitrary. Our presentation is as follows: in 



Sec. IIA we recall the basic strategy of the stiffness- 
matrix formalism considering constrained functionals of 
N collective coordinates and introducing the N ^~ 1 '> 
structure factor matri x elements S^ v (q) (with 1 < /a, v < 
N). Then in Sec. [I B we show how the local character of 
Sl[m] necessarily leads to two sets of algebraic relations 
relating the S^fa) and £^,,(0) at any three planes Z\, z 2 
and z 3 . This means that for arbitrary positions z\ < z 2 < 
z 3 the zeroth moments Gq(z^, z v ) satisfy two non-trivial 
identities. We illustrate this in Sec. Ill for the simplest 



possible scenario corresponding to fluid adsorbed at a 
single wall before turning our attention to the much more 
interesting problem of fluids confined in a parallel-plate 
geometry. For this case, if we consider examples in which 
the equilibrium one-body profile exhibits a reflection 
symmetry in the z — U plane, the algebra simplifies 
considerably. Even symmetry corresponds to profiles 
satisfying mo(z) = mo(L-z) and occurs if both confining 
walls are identical. For this case one of our relations 



2 



readily re-derives an exact result due to Henderson 
|L5| valid for fluids (interacting with quite arbitrary 
interatomic forces) confined between perfectly hard walls. 
On the other hand profiles satisfying mo(z) = —mo(L—z) 
correspond to odd symmetry and occur in Ising-like 
systems with competing surface fields. Such geometries 
have attracted considerable attention in recent years (see 
P,^| Jl^Jl7[| for example and references therein) due to 
presence of strong interfacial finite-size effects and novel 
symmetry breaking mechanisms. 

For both even and odd systems the correlation function 
relations severely restrict the form of the zeroth moment 
and we are able to derive a very elegant expression 
for Go(z\,Zz) in terms of Gq{z\, Z\), Gq(z2, z%), m' (z) 
and 7(L). Using this approach we rederive non-trivial 
scaling laws at the bulk critical point and when large 
scale interfacial fluctuations are present in the confining 
geometry. Sec. |v| concludes our discussion of the planar 
system where we consider the continuum (iV — ► oo) 
limit of the relations. We then turn our attention 
to non-planar geometries for which the stiffness-matrix 
formalism has not yet be fully developed. Nevertheless 
we show that for Landau-type models of inhomogeneous 
fluids in cylindrical and spherical symmetries it is 
possible to define analogous matrix elements S^ v which 
satisfy the same algebraic relations as those in the planar 
geometry. 

We conclude our article with a summary of our main 
results and make some remarks about other types of 
variational model. 



II. CORRELATION FUNCTION RELATIONS 

A. Constrained functionals 

Consider a fluid confined between two planar walls (of 
infinite transverse area A) situated in the planes z = 
and z — L. We consider the most general situation in 
which the walls are not identical and make no assumption 
about the preferential adsorption at each surface. Now 
consider an ordered arrangement of N planes (labelled by 
a = 1, . . . , N) corresponding to positions zi, Z2, ■ ■ ■ , Zjsr 
and denote the equilibrium magnetization 



m (z o 



(2.1) 



at each (see Fig. 1). We next introduce a constrained 
minimization of the functional fi[m(r)] which may be 
viewed as variational prescription for integrating out 
degrees of freedom except for those infinitesimally close 
to the planes. To this end suppose that fi[m] is minimized 
subject to the condition that the surfaces described 
by the position variables £ a (y) are contours of fixed 
magnetization 

m(4(y))=m* Va (2.2) 
Thus we define a constrained functional Jlif 




FIG. 1. Schematic illustration of the equilibrium profile 
mo(z) for a thin film geometry. The equilibrium positions of 
surfaces of fixed magnetization mf , and m$ are shown. 



H N [{£ a (y)};{z a }} = minO[m(r) 



(2.3) 



where the bar denotes the crossing criterion constraint 
( p.2[ ). By construction the equilibrium positi on of the 
surfaces of fixed magnetization is z a (see d2.l| )) and 
setting £ a (y) = z a in ( |2.3P recovers the thermodynamic 
potential Q. To make connection with the correlation 
functions we need to consider the properties of the 
constrained functional in the vicinity of the global 
minimum. To this end we suppose that for small 
translations and fluctuations we may write [[L2| 



H N [{£ a (y)};{z a }} = Jdy 

+w N ({e a y,{z a }) 



-£ M „({M;{> a })VVVC 



(2.4) 



where the E^j, constitute the elements of the stiffness- 
matrix X for this par ticular choice of \m ^} (or equiv- 



alently {z Q }, see ( |2.l| )). In writing (2.4) we have (for 
convenience) omitted to include terms related to the 
curvatures V 2 £ Q etc. which may be accounted for using 
an appropriate rigidity matrix K. Also in contrast to 
the analogous effective Hamiltonian expressions there 
is no need to specify the momentum cut-offs for each 
field. Obviously different choices of the {z a } will result 
in different binding potentials Wn and stiffness-matrices 
X. For this reason it is sometimes convenient to speak 
of a continuous set of functionals {H]y({£ a }; {z a })} 
whose elements are distinguished by the particular 
choice of {z a } or equivalcntly {m^}- Connection with 
the correlation function G(zi, Z2',q) is made with an 
appropriate finite dimensional matrix r epre sentation of 
the Ornstein-Zernike integral equation Jl.4| ), The 
analogue of the direct correlation function is an TV x TV 
matrix with elements (setting fc^T = 1) 



C^(yi 2 ;{z a }) = 



6 2 H N ({£ a };{z a }) 



^(yi)4(y 2 ) 



(2.5) 
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which we Fourier transform to exploit the translational 
invariance. The wave vector expansion of the matrix 
C(g; {z a }) is rather elegant 



C(g; {z a }) = C Q ({z a }) + g 2 £({z Q }; {z a }) 
+g 4 K({z Q };{z Q }) + ... 



(2.6) 



and is central to the utility of the stiffness-matrix 
formalism. The zeroth term Co is the matrix of 
curvatures 



C ({z a }) = 



r) 2 f) 2 
u ll °V. 

. d 2 



8 2 



di 



N 



W N ({£ a };{z a }) 



(2.7) 



where d 2 ^ = 



and is evaluated at equilibrium l a = 
z a . The correlation function G(z^,z v :q) then follows as 



G(z fl , z v \ q) = m' (z IJ ,)m' (z v )S fl , v (q; {z a }) 



(2. 



where the S^ v are the elements of the (symmetric) 
structure factor matrix satisfying 



S(g;{za})C(g; {z a }) 



(2.9) 



with I the identity matrix. These equations have 
previously been used to study mean-field correlation 
functions at wetting transitions flS]| . In that context the 
HN({£ a };{z a }), correspond to coupled effective Hamil- 
tonians and the minimization condition arises as a saddle 
point approximation to a functional integral. Here we 
emphasise that the formulation is an exact prescription 
for calculating correlation functions provided one begins 
with the appropriate variational functional. 



B. Local functionals 

To begin consider the zeroth moment Gq^z^zJ) for 
which we only need the binding potential. In order to 
derive the correlation function relations it is sufficient 
to consider the properties of 'three field' functionals 
H 3 \£\, £2, £3] ■ as working with N > 3 does not give 
any new results. The essential observation leading to 
the correlation function relations is that due to the local 



character of the grand potential functional (1.11) the 



binding potentials necessarily have a separable form (see 
Appendix A) 



W 3 ({£ a (y)};{z a }) = V 1 (£ 1 ;{z a }) 
V 2 (£ 2 - h; {z a }) + V 3 (£ 3 - £ 2 ; {z a }) 
+ V i (£ 3 -{z a }) 



(2.10) 



where the functions Vi depend on the particular choice 
of variational model and may be considered unknowns. 



From (2.7) it is a trivial exercise to calculate the matrix 
Co({ z a}) and note the following general properties 

Ci 3 (0;{z Q }) = (2.11) 
Ci 2 (0;{z Q }) + C 22 (0;{z Q }) + C 23 (0;{z Q }) = (2.12) 

In turn these impose conditions, via ( [2.9] ), on the 
structure factor matrices S(0;{z Q }) (hereafter we drop 
the explicit {z a } dependence). 



5i 2 (0) 5 13 (0) 
5 22 (0) 5 23 (0) 







and 



5x2(0) S 23 (0) 
5x3(0) S 33 (0) 



Sn(0) 5 12 (0) 
Si 3 (0) 5 23 (0) 



(2.13) 



5 n (0) 5 13 (0) 
5i3(0) 5 33 (0) 

(2.14) 



For q 7^ it is necessary to consider the properties 
of the stiffness-matrix and rigidity etc. appearing in 
the expansion ( |2.6[ ). Nevertheless because of the local 
character of f2[mj, variation of ^x(y), say, does not effect 
the constrained magnetization in the region z > £ 2 (y) 
(see Appendix A). Consequently the Si 3 element of the 
stiffness-matrix vanishes (as does i^i 3 the corresponding 
element in the rigidity matrix) leading to 



C 13 (q;{z a }) =0 Vg 



(2.15) 



which is clearly the generalization of ( 2.11 ) and leads to 
the first of our identities 



5i 2 (g)5 23 (g) = 5 22 (g)S* 13 (g) 



(2.16) 



consistent with ( [2.13 ) when g = 0. T he se cond identity 
is restricted to g = and is given in (2.14) above. This 
may be profitably rewritten as 



[5xx(0)- 5i2(0)] [5 33 (0) - 5 23 (0)] = 
[5i3(0) - Sia(0)] [5i3(0) - S 23 (0)} 



(2.17) 



If we were to consider an iV-field constrained functional 
Hn[{£oi}] {z a }] the S^i,(q) are related by the same 
relations applied to any ordered triplet of planes located 
at z a < zp < z-y, i. e. 

5a/3(g)5/3 7 (g) = Spp(q)S ay (q) 

and 

[S aa (0) - S a0 (O)} [5 77 (0) - 5^(0)] = 
[5 Q7 (0) - S aP (0)] [5 Q7 (0) - 5,3 7 (0)] 

Clearly ( 2.18| ) and ( 2.1 9| ) define an algebra satisfied by 

the jV ^~ 1) quantities S flv (q) (for each g). In terms of 
the pair correlation function G(z^, z„; g) they read 

G(zi, z 2 ;q)G(z 2 , z 3 ; q) = G(z 2 , z 2 ;q)G(z ll z 3 ] q) (2.20) 



(2.18) 



(2.19) 
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and 



with 



m' (z 2 )Go(z 1 ,z 1 )-m' (z 1 )Go{z 1 ,z 2 ) _ 
m' (z2)Go(zi .Z3) -m' (z3)Ga(zi ,zi) ~~ 

m' (z 2 )Go{z 1 ,Z3)-m' (z 1 )Go(z2,z 3 ) 
m' (z2)Ga(z3,z 3 )-m' a (z3)Ga(z2,z 3 ) 



(2.21) 



for all < z\ < z 2 < z 3 < L. We emphasise that these 
relations are valid for all local variational models of the 



form (1.11). As we shall see, taken together they restrict 
the form of the zeroth moment Gq(zi,z 2 ) and allow us 
to derive elegant expressions relating Go to the force of 
solvation f s {L). Before we do this we consider the case 
of flu id ads orptio n at a single wall for which the algebra 
( [2.18 ) and ( 2.19 ) has a trivial solution. 



III. A SIMPLE CASE: FLUID ADSORPTION AT 
A SINGLE WALL 

We wish to show how for a s emi- infinite system the 
algebra conditions, (fy||) and (|2~19| ), with q = are 
met. To this end we only need to consider the properties 
of the binding potential Wn{{(.^})- For later purposes it 
is convenient to position the wall in the plane z = zq 
(fixed). Clearly the equilibrium profile mo(z; zq) is a 
function of z — zq only so that the partial derivatives 



satisfy d z mo(z; zq) 



-d Zo m(z; zq). Due to the local 



character of the grand potential function the TV-field 
binding potential may be written 

Wn({^}) = v x {h) + V 2 (£ 2 - ti) + v 3 (e 3 - l 2 ) 

+ ... + V N {t N -i N . l ) (3.1) 

The simplifying feature here (compared to the parallel- 
plate geometry) is that the final collective coordinate 
•^zv(y) only enters through one (unspecified) partial 
binding potential function Vn(x). The curvature matrix 
has the tridiagonal form 



-V," 



v 2 " 



v 3 " 



■vi' 



( v{' + y 2 " 

-v 2 " 





v ; 



where we have written V' 



\ 



-V" 

* N—l 



-V" 



v N 



— V" 
V" 



d 2 V a (x) 



equilibrium 



(3.2) 

(evaluated at 



The inverse matrix Sq has a 



remarkably simple block structure 



So 





Sn 


Sn 


Su 




Su 


Su 


s 22 


s 22 


s 22 






S u 


s 22 


S33 


S33 






Su 


s 22 


S33 










s 22 


S33 






Sn-2 : N-2 










Sn-i,n-i 


Sn-i,n-i 


\Su 








Sjv-1,N-1 


Snn 



(3.3) 



Sfi/j, — <s^(o) 



1 



a=l <* 



(3.4) 



The point to notice here is that 

V(0) = 5 W (0) -iv>ii (3.5) 

which is i ndeed a solution of the algebraic conditions 
we are led to the conclusion that within 
all local variational models of fluid adsorption at a single 
wall the pair correlation function is of the form 

G (zi,z 2 ) = d z m (zu z Q )d z m (z2; z )F (min[zi, z 2 \) 

(3.6) 

or equivalently 
G {zi,z 2 ) = d Za m (zi; z )d Zo m (z 2 ; z )F (min[zi, z 2 }) 

(3.7) 

where F(x) is an unknown function. 

This prediction is supported by the explicit Landau 
theory result for a semi-i nfinit e system (i. e. if we assume 
that £W is given by ( |l,12|) ) which is known to be 

G (zi,z 2 ) = m' (zi)m' Q (z 2 ) I a + l 



with 




m (0)Km (0)-m '(Q)] (3.9) 
d z mo(z; zq) and 



where we have abbreviated Tn' (z) 

Equation ( |3.5| ) is also consistent with an exact result 
due to Henderson and van Swol p0[ who have derived an 
exact statistical mechanical sum rule for fluid adsorption 
at a pure hard wall i. c. an external potential 



Vhw(z) = 



00 




for z < 
for z > 



(3.10) 



which simply confines the fluid to the half plane z > 0. 
For this case it is possible to show that 

G o (0, z) = p' (z) for hard walls (3.11) 

with po(z) mo(z)) the equilibrium local number 
density. This is in agreement with ([Tf]) for the special 
case z^ — equivalent to 

G o (0,z) cx m' (z) (3.12) 

with an unknown constant of proportionality which 
presumably depends on the choice of surface interaction 
term <pi(m). 



While the semi-infinite solution (3.5) is not particularly 
interesting it does represent the correct L — > 00 limit of 
the algebra pertinent to fluids confined in parallel-plate 
geometries. The analysis for this case is much richer and 
is discussed in the next section. 
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IV. PARALLEL PLATE GEOMETRIES 

A. Preliminary remarks 

There are t wo ca ses w hen the correlation function 
relations ( 2.1S| ) and ( |2.19 ) — with q = 0, simplify further 
and yield elegant expressions for Gq{z\,Z2). These 
occur in parallel-plate geometries in which the one- 
body profile Too (z) exhibits a simple reflection symmetry. 
For example if both walls are identical it necessarily 
follows that too {z) is an even function about the plane 
of symmetry z = |-. These systems have traditionally 
attracted the most attention in the literature where the 
shift of the bulk critical point and first-order phase 
boundary (capillary condensation) are of interest. The 
finite-size phase diagram and representative profiles for a 
typical parallel-plate geometry are shown schematically 
in Fig. 2. However more recently, examples in which 



fTc 
T C (L) 



m (z) 



T 

A 




h 



h c (L) 

FIG. 2. Schematic phase diagram of a thin film Ising-like 
magnet with positive surface fields hi = h.2- Representative 
magnetization profiles, shown inset have an even reflection 
symmetry. The locus of first-order phase transitions (often 
referred to as capillary condensation) is shifted away from 
the zero bulk field (h — 0) line. The shifted capillary critical 
point occurs at Tc{L) and hc{L). 



the profile exhibits an odd reflection symmetry have also 
drawn considerable interest. Such parities arise in Ising- 
like systems confined by walls which exert surface fields 
hi and hi of equal magnitude but opposite sign on the 
spins in the z = and z = L planes respectively. 
The nature of the phase coexistence and criticality in 
the system is entirely different to that occurring for 
the case of identical walls. In particular the critical 
temperature Tq(L) for finite L is determined by length 
scales associated with wetting |jj and is restricted to lie 
close to the wetting temperature Tw of the semi-infinite 
system (see Fig. 3). For temperatures T > Tc(L) and 



One phase 
region 



Two phase 
region 



Tc 





m (z) 









L 















T C (L) 



h 





FIG. 3. Schematic phase diagram of thin film Ising model 
with opposite surface fields hi — —h^. For temperatures 
T > Tc(L) (and zero bulk field) the magnetization profiles 
have an odd reflection symmetry. In the temperature window 
Tc > T > Tc{L) the profile resembles an interface located 
near the middle of the thin film and characterized by an 
extremely large interfacial correlation length . The location 
of the capillary critical point is determined by lengthscales 
pertinent to wetting and so Tc(L) is less than but close to 
the critical wetting temperature Tw- 



h — the magnetization profile has the odd symmetry 
specified above. Of particular interest is the temperature 
window Tc > T > Tw corresponding to the regime 
where the finite-size effects have suppressed bulk phase 
coexistence. The profile Too {z) for these temperatures 
resembles an up spin-down spin interface situated in the 
centre of the system. This interface is very weakly pinned 
by the confining walls and wanders almost freely in the 
finite-size geometry. The fluctuations associated with 
this wandering are extremely large leading to universal 
scaling behaviour for sufficiently low dimension p|. We 
shall return to this in Sec. IV C where we use the relations 



derived in Sec. IV B to calculate the singular contribution 
to the finite-size free-energy j(L) (for both even and odd 
systems) . 

To begin however, we make some prelimary remarks 
which will be useful for later purposes and allow us to 
make contact with other work. To establish our notation 
we note that the magnetization profile satisfies 



TOq(z) = ±m (L — z) 



(4.1) 



for even (— ) and odd (+) systems. Moreover we can also 
impose the additional symmetry requirement 

G (z l ,z 1 ) = G (L-z 1 ,L-z 1 ) (4.2) 

or equivalcntly 

Sn(0) = S 33 (0) {ovz 3 =L- Zl (4.3) 



G 



Substitution into ( |2.19D yields the 'sum' rule 



B. Connection with the free-energy 



Si2(0) + S 23 (0) = Sii(0) + Si 3 (0) z 3 = L- Zl (4.4) 



which compares to the 'product' rule ( 2.1 8| ) 

5i 2 (0)5 23 (0) = S 22 (0)S 13 (0) (4.5) 



We can therefore anticipate that S\ 2 (0) and S* 23 (0) must 
correspond to the roots of a single quadratic equation 
(see S ec. | IVB| below). In terms of the zeroth moment 
itself (}4.4[) implies 



G (z 1} z 2 ) ± G {z 2 , L - zi) _ m' (z 2 ) 
Gq{zx,z\) ± G (z 1 ,L- zi) m' (zi) 



(4.6) 



for even (— ) and odd (+) systems and z\ < z 2 < L— z\. 
Further insight follows if we set Z\ = to find 



G (0,z)±G o (z,L) = cm' {z) 



(4.7) 



where c is an unknown constant (which may depend on 
T,L,h,...). Note that the equilibrium magnetization 
profile is also a function of these variables. 

It should be emphasized at this stage that this elegant 
relation follows as a necessary consequence of th e local 
nature of the underlying variational model (1.11). It is 
therefore encouraging to note that again it is consistent 
with an exact result due to Henderson |fl5|| who has 
considered statistical mechanical sum rules for a fluid 
with arbitrary intcrmolecular forces confined between two 
identical hard walls (i. e. even symmetry). Henderson 
supposed that the external potential could be written as 
a sum of two semi-infinite wall contributions 



V e ^(z) = V 00 {z) + V 00 (L-z) 



(4.8) 



and derived a number of exact relations between integrals 
over Gq(zi, z 2 ) and one-body/thermodynamic quantities. 
For the specific case of hard walls Voo — Vrw these 
simplify and in particular yield 

G (0, z) - G (z, L) = p' Q (z) for hard walls (4.9) 



with po(z) rno(z)) the equilibrium (local) number 
dens ity. Clearly this is in agreement with our prediction 
(4.7) and identifies the unknown constant c a s th e 
universal value chw = 1> for hard walls (3.10). 
Henderson [fL5| also derives a relation between the pair 



correlation function and the free-energy for this system 
(recall that we have set k^T = 1) 



for hard walls 



(4.10) 



In the next section we show how this (and more) can be 
derived from the correlation function algebra for more 
general wall potentials. 



For q — there are two identities satisfied by the 
S^i^O) at any three planes z\ < z 2 < z 3 valid for 
arbitrary (short ranged) wall interactions cf>i and </> 2 . 
To develop the theory further beyond the elementary 
remarks made above, it is necessary to inquire what 
restrictions these conditions imp ose o n the structure of 
5^1/(0). Consider for example ( 2.18 ), with q = 0. In 
order that this is satisfied for arbitrary choices of z\ , z 2 
and Z3 it follows that the structure factor must have the 
form of an ordered product 



for z p < z v 



(4.11) 



regardless of whether the profile mo(z) exhibits a 
reflection symmetry or not. Here y~{z) and y + (z) are 
unknown functions, the properties of which we need to 
determine. Note that the ordering condition z M < z v 
is of some importance here and is less restrictive than 
the assertion that ^^(0) is a separable function of z M 
and z v . While it is clear that similar remarks also apply 
for q ^ 0, we do not have a second relation which yields 
valuable information about y~{z) and y + (z). To proceed 
we substitute ( [4.1 fl) into ( |2.19| ) to find 



y + (zi)y [z 3 )-y + (zi)y (z 2 )-y + (z 2 )y (z 3 ) = 
y+(z 3 )y~( Zl ) - y-(z 2 )y+(z 3 ) - y+(z 2 )y-( Zl ) (4.12) 

which can be differentiated with respect to z 2 yielding 
(for all z\ and z 3 ) 



(4.13) 



,2/ O3) -y (zi) 
y-(z 2 )' = f3y + (z 2 )' 



for constant /?. This relation is valid for all z 2 and can 
be integrated to show that y~(z) and y + (z) are linearly 
related 



y + {z) 



(4.14) 



with a the constant of integration. By using (4.11) and 
( fl.l4| ), we can solve for y _ (z M ), say, in terms of S^^O) 



-V(O). We find 



and then use (4.11) again to yield an expression for 



KS fiV (0)= \1±J1-KS, 



(0) j (l ± y/1 - KS vv (ofj 
(4.15) 



where K = — ^ is a single undetermined constant. The 
± signs do not correspond to the odd/even reflection 
symmetries mentioned in Sec. IV A since we have not 



yet specialised to these systems. Instead they refer to 
the positions of z^ and z v relative to the maximum of 
Gq(z, z), as we shall see below. 
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Before we turn to odd/even systems, for which we can 
explicitly calcula te th e constant K, we note that the 
general solution (4.15) is consistent with remarks made 
earlier concerning the semi- infinite limit. The alge bra for 
this case corresponds to the K — ► limit of (4.15), with 
the appropriate choice of signs, 



variables at any two planes z^ and z v {> z^) are related 
by 



= (1 ± y/l - 5~) (1 ± y/1 - o vv ) (4.22) 

where the ± signs must be chosen appropriately (see Fig. 
4). It is a straight forward exercise to check that this 



lim 



o 5^1,(0) 



lim ] , 



1 - KS^(0) 



y/1 + KS VV {Q)) 



V(o) 



(4.16) 



as indicated in (3.5). 

Hereafter we speci alise to odd (+) and even (— ) 
syste ms. Rearranging ( 4.15| ) and using the results of Sec. 
IV A it is straight forward to derive 



K± 



±AG Q (0,L)m' Q (0) 2 
(G o (0,0)±G o (0,L)) 2 



(4.17) 



However it is also possible to relate K± to the excess 
free-energy "f(L). We simply quote the result and refer 
the interested reader to appendix B for the details 



K± = 4 



d 2 l(L) 
dL 2 



(4.18) 



For odd symmetric systems it is also straight forward to 
derive 



HH 2 



2 I 



G ( 2", g) 



(4.19) 



which is particularly useful. Note that the right hand 
side of this relation can be identified as the reciprocal 
of the maximu m va lue of S uu (0 ) as a function of 



In fact, using ( 4.11 ) and ( 4.14j) to write S MM (0), and 



assuming analyticity at all points (true for odd systems 
only), simple differentiation with respect to z M shows the 
maximum value to be 5^(0) = — fg, which is just j^. 

Thus we have established the following relationship 
between the correlation function and excess free-energy 



d 2 1 
dL 2 



± 



G o (0,L)m' o (0) 2 
[G o (0,0)±G o (0,L)] 2 



(4.20) 



which is once again (remem bering (4.7 )) consistent with 
Henderson's exact r esult s (4~9|,4.10) for fluids confined 
between hard walls ( 3.1 Cl| ) . 

In this way we are led to a very elegant universal equa- 
tion relating suitably scaled moments of the correlation 
function. Specifically, define the dimensionless quantities 



4 d 2 7 Gq{z, 1 ,z v ) 
k B T dL 2 m! Q (zi l )'m! Q {z v ) 



(4.21) 



where for completeness we have reinstated the Boltzmann 
factor. For both odd and even systems the value of these 




L/2 L 

FIG. 4. Choice of signs appearing in the first and second 
parenthesis, respectively in equation (4.24). 



result for a^y is consistent with the algebraic relations 
CT12O23 = er22Ci3 (4-23) 

((Til - 0-12X033 _ °23) = (013 - Ol 2 )(o-l3 - cr 23 ) (4.24) 

for z\ < Z2 < Z3. 

Making note of the sign of K± we see that even and 
odd systems are distinguished by the ranges of values of 
the nil variables. In particular 

< o^v < 1 for odd systems (4-25) 

with <7 MM = 1 for z,j, = j. On the other hand in even 
systems where the force of solvation f s (L) between the 
plates is attractive we have 



07,^ < for even systems (4.26) 

and is not bounded from below. Note that <7 W — > — 00 as 
z^ — > i due to the turning point in the even equilibrium 
magnetization profile (sec Fig. 2). For both even and odd 
systems the variables ct m „ and ct ma1 approach zero when 
L — > 00 (for fixed temperature T) consistent with the 
semi-infinite limit K± — > mentioned earlier. 



Equations (4.17-4.22) are the main results of this paper 
and are obeyed by all local variational models. In the next 
section we shall demonstrate how they may be used to 
derive non-trivial scaling expressions for the fluctuation 
contribution to the free-energy j(L). 



C. Scaling of the finite-size free-energy 

There are three different scenarios involving fluids 
confined in parallel-plate geometries for which the free- 
energy j(L) may be expected to exhibit singular, 
non-analytic behaviour arising from fluctuation related 
effects. These arc considered separately below 
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1. Odd and even systems at bulk criticality 

Exactly at the bulk critical point T = Tc and h = the 
correlation length in the finite-size system is only limited 
by the width L. This implies long-ranged (i. e. algebraic) 
behaviour in the finite-size free-energy j(L) which may 
be likened to the Casimir effect in quantum field theory 
pl[ . The required non-classical (i. e. non-mean-field) 
scaling of 'y(L) may be derived from our equations in 
a number of ways. Perhaps the most elegant is to assert 
that due to the universal character (i. e. temperature and 
bulk field inde pendence) of the scaled correlation function 
relation ( 4.22 ) the variables cr M „ must exhibit finite-size 
scaling of the form 



(4.27) 



where the scaling functions ^^(x, y) are universal |2j ] 
for all odd (+) and even (— ) systems (restricting our 
attention to surface fields hi = ±/i2(> 0) and positive 
surface enhancement). This hypothesis may be compared 
with the definitions (4.21) involving the free-energy. 
Standard finite-size arguments dictate that sufficiently 
far away from the wall the correlation function and 
magnetization profile behave (exactly at criticality) as 



and 



G (zi, z 2 ) 



m {z) 



(!) 



(4.28) 



(4.29) 



where g and A are appropriate scaling functions, and 7, v 
and P are s tand ar d bul k critical exponents. Combining 
these with fl4.2l|), fl4.27|) yields 



7 (L) ~ L 1 - 



at T = T c 



(4.30) 



where we have used the Rushbrooke relation 2 — a = 
2(3 + 7. Thus for d > 4 we expect j(L) ~ L~ 3 while for 
d < 4 we find 'y(L) ~ L~ on invoking hyperscaling. 
This is precisely the predicted scaling behaviour of the 
surface excess free-energy |l|,p||2ll| and serves to ill ustrate 
the generality of the relations derived in Sec. [V B . While 



they may be of restricted validity, applicable only to local 
variational model, they are not necessarily mean-field-like 
in character. 



One may also derive ( 4.30 ) using ( 4.2C ) which implies 



^«G o (0,L) 



as L 



(4.31) 



at all temperatures. Using the 'product' relation one has 
G o (0,L) = ^TO (4.32) 



G o (j> f) 



and it is then straight forward to ca lcula te the L 
dependence of the right hand side of ( 4.31 ) at bulk 
criticality. The denominator follows from (4.28) yielding 
Go (-§,■§) ~ L 77 " 1 . In principle the scaling of Go(0, -§) 
could be calculated fr om t he short distance expansion 
(valid for 4^ <C 1) of ( 4.28 ). However it is easier to use 
the semi-infinite result Go(0,z) cx rn' (z) quoted in Sec. 



III. Combining these with the finite-size scaling of the 



profile ( 4.29[ ) it is natural to suppose 



G o (0,^)~L 
which recovers j(L) 



at T = T c 



(4.33) 



-+1 



directly without 



appealing to the finite-size scaling of <t m „ . 



2. Odd systems with large interfacial fluctuations 

A second example in which the free-energy exhibits 
long ranged behaviour occurs in odd systems in the 
temperature window Tc > T > Tw (often referred to 
as the soft mode phase [|§l]4| , 0, |i^ , |r^ . For this case it is 



most convenient to use ( ]4.18[) and ( 4.19 ) 



dL 2 



G (, , TT J 



(4.34) 



and recognise that the right hand side is proportional 
to S,^ 2 , where £|| is the transverse correlation length 
characterising fluctuations in the position of the up spin- 
down spin interface, which wanders 'freely' between the 
two walls. If we recall this correlation length shows finite- 
size scaling such that (J] 



(4.35) 



with C the roughness exponent we immediately find 



7 (L) ~ L 2 * 1 "*) 



(4.36) 



consistent 
interfaces 



with the fluctuation theory for confined 
22j . For purely thermal fluctuations for which 



C = for d < 3 this reduces to 



7(L) - L— 



for d < 3 



(4.37) 



familiar from the theory of wetting |2^]. The same 
expression also follows using arguments similar to those 
given for the case of T — Tc and serve to further illustrate 



the utility of the results (|4 . 1 7| - |4 . 2 2| ) 



3. Near the finite-size critical point 

Finally we make some remarks on a third possibility 
which requires further research. For both odd and even 
geometries it is natural to expect singularities to emerge 







as we approach the finite-size critical point occurring at 
Tq(L) and h = hc(L). Recall that the locations of 
these critical points is very different for odd and even 
symmetries (see Fig. 2 and Fig. 3). Nevertheless the 
phase transition occurring at Tq{L) and h = hc{L) for 
both these systems is conjectured to belong to the same 
d — 1 dimensional bulk Ising universality class [p|,pl|. 
However an inspection of the well developed mean- 
field theories shows that the all important free-energy 
derivative ^p- has different singularities for odd and even 

,2 

symmetry In particular diverges at the finite-size 
critical point for an even system jl9| but vanishes (as 
Tc -> T C (L) + ) for an odd symmetry j|. This would 
seem to imply a subtle difference between the behaviour 
of the correlation functions near the respective critical 
points (by virtue of (4.18) and (4.22)) — an observation 
that is perhaps surprising given, as mentioned above, 
that the universality class of the phase transition is 
anticipated to be the same for each geometry. 



V. ALTERNATIVE APPROACHES AND SOME 
GENERALIZATIONS 

To end our discussion of planar inhomogeneous fluids 
we present an alte rnative d erivation of the correlation 
function relations (2.18 2.19) for Landau- type models of 
the form 



Q[m(r)] = J dy J dz ji(Vro) 2 + <f>(m) + S(z)<j>i{m) 

+6(z-L)<h(m)} (5.1) 

by p artial solution of the Ornstein-Zernike equation 
(1.7) — a complete solution not being required. Here 
we emphasis the properties of a kernel or propagator- 
like function which also emerges by considering the 
cont inuum limit of the algebraic relation fl2.18| ) (see Sec. 
VB). Finally we demonstrate that the same algebraic 



relations amongst the (suitably redefined) Sa V ar e also 
obeyed for non-planar inhomogeneous fluids exhibiting a 
cylindrical or spherical symmetry. 



A. The Kernel Function for the Ornstein-Zernike 
equation 



For the Landau model (5.1) the Ornstein-Zernike 
equation reduces to 



£(z 2 ) + q 2 G(zi,Z2',q) = S(z 2 - z\) 



where the second order linear operator C is 



C(z) = -— + q>"(mo(z)) 



(5.2) 



(5.3) 



The same operator appears in the differential equation 
for the profile ttiq{z) which may be written 



C(z)m' (z) = 



(5.4) 



and will be required later. To proceed we first define the 
function 



d 

M«i,22;q) = ^iogG(zi,z 2 ;q) 



(5.5) 



which from (5.2) satisfies the non- linear equation (for 
zi ^ z 2 ) 



^ = Am(z 2 )) + g 2 -/C 2 

dZ2 



(5.6) 



Assuming (5.6) has a family of solutions fc(z 2 ,q;a), 
parameterized by a, then 

}C(zuZ2;cO = { k S Z2 ^ ai \ f TZ2 >* (5.7) 
v 1 fc(z 2 ,q;a 2 ) for z 2 < z\ v ; 



and the 5 function in (5.2) implies the boundary 
condition 



fc(zi,q;ai) - fc(zi,q; a 2 ) 



1 



G(zi,zi;q) 



(5.8) 



Without loss of generality this can be considered an 
equation for a 2 and implies that a 2 = a 2 (zi,q). 
Consequently, 



1C(zi, Z2] q) = 1C(Z2] q) for z\ < z 2 



(5.9) 



As we shall see the existence of K. is central to the deriva- 
tion of the correlation fun ction relations. Integrating 
Q5.5D and remembering ( |5.9D we find (for Z\ < z 2 ) 

G(z 1 ,L;q L ) - G(z 1; z 2 ;q) = / dz'JC(z'; q)G(z' , Z\\ q) 

J Z 2 

(5.10) 

This becomes in real space (via the convolution theorem) 

G(r 1 ,r 2 ) = G(r 1 ,(y 2 ,L)) + J dr'G(r u r')K-(r', r 2 ) 

(5.11) 

where z\ < z 2 and K~ can be thought of as an advanced 
propagator 



K-(r',r 2 ) 



{z - z 2 )K(z';y 2 - y') 



(5.12) 



Here K(z; y) is the Fourier transform of the kernel 
/C(z;q). In this way it is possible to interprete the 
correlation function at positions ri and r 2 (say) in a 
standard way i. e. making analogy between position z and 
time, G(ri,r') at early 'time' z = z\ can be propogated 
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to later time z = z 2 making use of information from the 
'future' [z! > Z2) only. 

The 'product' rule foll ows d irec tly from the integral 
relation (found from (5.5) and (5.9)) 



G(zi ,b;q) — G(zi , a; q) exp 



dzfC(z; q) 



(5.13) 



valid for all b > a > z\. Setting z\ = a and choosing c 
satisfying a < c < b we can write 



G(a, b; q) = G(a, a; q) exp 
= G(a, a; q) exp 



x exp 



G(a, a; q 



dzJC(z; q) 

dzK.(z; q) 

J dzfC(z;q 

G(a,c;q) G{b, c; q) 



G(a, a; q) G(c, c; q) 



(5.14) 



where we have used (5. 13 ^extensively. The final equation 
is nothing more than ( [2.1£ ) . 

The second correlation fun ctio n relation can be found 
from the explicit solution of (|5.6|) when q = 



/C 



dz 



[io g y (2)] 



(5.15) 



, = ^( 

one finds (for z 2 , z 3 > 2j ) 



where = ^"(rao(2))^( z )- From the definition ( |5.5| ) 



■G (z 1 ,z 2 ) = Y(z 2 )X(z 1 ) 



(5.16) 



for some function X(z). Enforcing the Ornstein-Zernike 
conditions 

C(zi)G (zt, z 2 ) = C(z 2 )Gq(zi,z 2 ) = S(zi - z 2 ) (5.17) 

we see that the functions X and Y are related by 

X(z) = aY{z) +(3m' {z) (5.18) 

for constant a and /?, due to the linear nature of the 
operator £(z) and (|5.4l) . This is precisely equivalent to 
the linear relation Q4.14| ) which we deduced as a necessary 
conse quence of the second correlation function identity 



B. The continuum limit 



The propagator approach also emerges from the 
more general stiffn ess- m atrix formalism by taking the 
continuum limit of (2.18|). Consider two boundary planes 



located at z = a and z — b (with a < b) and divide the 
region between them by n other planes separated by a 
constant distance h. Combining 'product' rules ( 2.18| ) 
defined for each consecutive set of two planes and the 
boundary at z — b one finds 



G(a,6;q) 



G(a, z 1 ; q)G(zi , z 2 ; q) . . . G{z n , b; q) 
G(zi,zi, q)G(z 2 , z 2 ; q) . . . G(z n , z n ; q) 

(5.19) 



Setting zq = a and z n+ i — b, ( 5.1 9| ) can be rewritten as 

G(a,b; q) _ ni=o z i + h \ q) 

G(a,o;q)~ ^=0 Zf, q) 

= exp[Sf =0 {log G(z,, + ft; q) 
-logG(^,^;q)}] 

logG(zj,Zj + ft;q) 



exp 



logG(z»,z»;q) 
ft 



(5.20) 



Taking the continuum limit ft — > 0, n — > 00 we find that 
G(a,6;q) = [/(a, 6;q)G(a,a;q) for a < 6 (5.21) 
where 

r 6 



C/(a, 6; q) = exp 



dzlC(z; q) 



(5.22) 



and /C(z; q) is some unknown function. For an infinitesi- 
mal displacement dz 2 , 



U(zi,z 2 + dz 2 ; q) = exp 



z 2 +dz 2 



dz fC(z; q) 



exp[fC(z 2 ;q)dz 2 ]U(zi,z 2 ; q) 
[l+IC(z 2 ;q)dz 2 ]U(z 1 ,z 2l q) (5.23) 



which implies 

U(z\, z 2 + dz 2 ) - U(zi,z 2 ) 
dz 2 

d 



= ]C(z 2 ;q)U(zi,z 2 ;q L ) 



— U{zi,z 2 ;q) =K,(z 2 ;q)U(zi,z 2 ;q) (5.24) 
oz 2 



taking the limit dz 2 — ► 0. It is intriguing that (5.24) 
is very similar to the Schrodinger equation for the time 
evolution operator, with the Hamiltonian time dependent 
but self comm uting at different times U~~ 



Multiplying ( 5.24 ) by G(zi, Zi; q) and then using ( 5.21 ) 



d 

dz 2 



U(zi,z 2 ;q)G(zi,zi;q) = )C(z 2 ;q)U(zi,z 2 ;q) 

xG(zi,zi;q) 

d 

=> —G(z 1 ,z 2 ;q) = K(z 2 ; q)G(z u z 2 ; q) (5.25) 
oz 2 
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which is equivalent to ( |5~5| ) with (5.9), and /C(z;q) is, of 
course, just the kernel function. 

Thus we have shown that ( [2.18 ) implies via its 
continuum limit the existence of a kernel function 
K,(z; q). In turn the kernel function can be understood as 
being an advanced propagator for correlation functions in 
real space. 



C. Non-planar geometries 

The properties of inhomogeneous fluids confined in 
non-planar geometries has also attracted attention in 
recent years. Perhaps the most important of these is the 
ease of fluid adsorption in cylindrical systems (as 
idealised models of porous materials). Despite continued 
interest we are not aware of any discussion of correlation 
function structure in such systems and to complete our 
article we make some remarks which follows from those 



mentioned in Sec. V A 



Consider then a fluid confined in an infinitely long 
cylinder (—00 < z < 00) of radius R (and connected to 
an external reservoir of particles). The grand-potential 
density functional is taken to be the Landau-type model 

/oo />27r rll r y 

dz I d(f> rdr < — (Vm) 2 + cf>(m) 
-00 Jo Jo 12 

+5(r - R)M m )} (5-26) 



analogous to (5.1). The connected correlation function 
G(ri,r 2 ) only depends on the relative angle 4> 2 — 4>i = 
021 and azimuthal distance 22 — z\ = z 2 \ between the 
particles as well as the radial distances r\ and r 2 . To 
exp loit this we first define a 'zeroth' moment analogous 
to (O) 



G (ri,r 2 ) 



dz 2 i 



271 



*2iG(ri,r 2 ) (5.27) 



The Ornstein-Zernike equation is 

[V 2 + <j>"(m(r 2 ;R))] G(r 2 ,n) = S(r 2 - n) (5.28) 

where m(r; R) is the equilibrium profile satisfying the 
Euler-Lagrange equation 

~ (r^j m(r; R) = <t>'(m(r; R)) (5.29) 

together with a boundary condition at r = R. The 
differential equation for G(ri, r 2 ) can be integrated with 
respect to z 2 \ and <p 2 i to yield (for r% 7^ r 2 ) 



1 


d 




r 2 


dr 2 





■4>"[m{r 2 ;R)) 



G (n,r 2 )=0 

(5.30) 



which we need to solve. The point to notice here is that 
the partial derivative 
differential equation 



the partial derivative ^g(r;i?) satisfies the same linear 



1 


d 




r 2 


dr 2 





+ cf ) "(m{r 2 ;R)) 



dm 
OR 



(r 2 ;R) = 
(5.31) 



As a consequence it is natural to define the 'structure' 
factor matrix elements as 



5^(0) 



OR ( r /ii R) gji (Tin R) 



(5.32) 



which, we anticip ate, s atisfy the same algebraic condi- 
tions ( ^.18 ) and ( 2.19 ) as the planar system. This is 
borne out by the explicit solution 



Otti dm 
Go(r 1 ,r 2 ) = M (r 1 ; J R) — 



(r 2 ;R)(a 



drr 



( dm 
VdR 



(5.33) 



where a is determined by a simple boundary condition. 
In terms of the structure factors we note that 



(5.34) 



-WO) = -WO) V ?V > r. 



which is analogous to the semi-infinite solution (|3.5| ) 
for planar systems. Indeed for fluid adsorption on the 
outside of a cylinder the Landau theory result is 



G (n,r 2 ) 



dm dm 
dR {n]R) dR 

min(ri,r 2 ) 



(r 2 ;R)(a' 

drr^ 1 1 — — 
\8R 



(5.35) 



which is precisely of the form (^ 

It is a simple matter to repeat the analysis for fluid 
adsorption in spherically symmetric systems (coordinates 
r, 9, <ft) where the appropriate zeroth moment is 

/>7T /■27T 

G (ri,ra)= / sm(6 12 )d6 12 # 2 iG(r!,r 2 ) (5.36) 
Jo Jo 

where 9 2 — 9\ = 9\ 2 . In fact, if we restrict our attention 
to adsorption at a single wall (with z <-> r and Zq R), 
or outside cylinders and spheres the general solution can 
be written 

„ , , dm . dm . ^, / 

G {n , r 2 ) = — (ri ; R) — {r 2 ; R) (^constant 

min(ri f ra) / q 

- -d' l ° 



^ I ~dR I I (5 ' 37) 



where d! = 0,1,2 for walls, cylinders and spheres 
respectively. We can interpret d! as the dimension of 
the boundary surface (two for three dimensional systems) 
minus the number of dimensions in which the systems 
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is unbounded. F or all these systems the structure 
factors satisfy (3.5). For the more general case of thin 
film geometries and adsorption between two concentric 
cylinders or spheres the structure factors can be s hown 
|^| to satisfy the general algebraic conditions (2. IS) and 



VI. CONCLUSIONS 

To complete our article we summarise our main results 
and make some remarks about possible future work. 

(i) The stiffness-matrix formalism for calculating cor- 
relation functions G{z ll ,z u \q) in planar inhomoge- 
neous fluids modelled by local variational models 
naturally leads to two identities relating this 
function at three arbitrary positions z\ < z 2 < 
23. Taken together they restrict the form of the 
zeroth moment Gq{z^,z v ) and allow us to define 
dimensionless scaled variables 



4 rf 2 7 Gojz^Zy) 
k B T dL 2 m'^z^m'^Zv) 



(6.1) 



for thin film geometries exhibiting even or odd 
reflection symmetries. The cross term a\ 2 is solved 
for explicitly as a universal function of an and a 22 . 

(ii) For even systems and particular choices of z\ , z 2 
and Z3 our results are consistent with previously 
derived exact expressions for fluid adsorption at 
purely repulsive hard walls |H| with arbitrary fluid- 
fluid interactions. 

(iii) Using our results we have been able to rederive 
non-trivial scaling expressions for the long-ranged 
(power law) finite-size contribution to the excess 
free-energy -f(L) 

(a) at the bulk critical point 

(b) in the soft mode phase for a geometry with 
competing surface fields (odd systems) 

There is a hint in our analysis that the correlation 
function may behave somewhat differently near the 
capillary critical point of odd and even symmetric 
systems. 

(iv) While the correlation function identities may also 
be derived directly, for Landau-type models, by 
explicit solution of the Ornstein-Zernike equation, 
it is unlikely that the general character of the 
relations would have been spotted using this 
method. This serves to illustrate the utility 
of the stiffness-matrix formalism which has been 
previously used to derive the stiffness-matrix free- 
energy relation in the theory of wetting |ll| . One 
may add here that while most of our analysis has 
been restricted to q — one may also discuss 



the structure of higher moments of the correlation 
function. For example, consider the position 
dependent transverse correlation length £||(z M , z v ), 
defined via the asymptotic expansion (for q — > 0) 



G{z^,z v ;q) = G (z^,,z u )[l - v)g 2 

+0(q 4 )] 



(6.2) 



Then, from the 'product' identity ( 2.1S| ) we im- 
mediately derive the elegant correlation length 
relation 

Z\\(zi,z 2 ) + ffi (z 2 , z 3 ) = £y (z 2 ,z 2 ) + Cf(«i s «a) 

(6.3) 

for z\ < Z2 < 23. Using this identity and other 
stiffness-matrix relations it is possible to show |26J] 
that the three correlation lengths £m (j, -§)' 
and £ 2 (0, L) all diverge in precisely the same 
manner as L — > 00 at T = Tq (even and odd 
systems) and Tq > T > T\y (odd systems) in zero 
bulk field. 



(v) Finally, we believe that it would be worthwhile 
developing the stiffness-matrix formalism for other 
variational functionals, for example, one could con- 
sider local models of two or more density variables 
such as arise in systems with iV-component vector 
order parameters. Within the same framework 
it would be possible to study local entropy-type 
functionals Q[m(r), e(r)] of the magnetization and 
the energy density which have been forwarded 
recently as better candidates for modelling critical 
behaviour po| . A preliminary analysis of these 
models reveals analogues of the 'product' rule 



(2.18) although further work is required to find 
whether the identity (2.19) also generalizes. 
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APPENDIX A: A SEPARABLE BINDING 
POTENTIAL 

Let us denote mi (z; {£ a }) as the planar profile which 
minimizes £l[m] subject to the crossing constraints ([2.2]). 
The binding potential is determined by the relation 
(dropping the implicit {z a } dependence) 
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W N ({£ a }) = Mmi N) (0; {*«})) + h(mi N) (L; {*«})) 
+ J L dzC® (mi N \z; {£ a }), d z mi N \z; {£,})) 

neglecting constant terms independent of the collective 
coordinates. This is profitably rewritten 

W N ({£ a }) = Mmi N) (0;{£ a })) 



+ Mmi N) {L- {M)) + E«=o ft* 1 dz 
(mi N \z; {£ a }), d z mi N \z; {£ a }) 



(Al) 



where we have defined £o = and £n+i = L. The profile 
itself satisfies the second order Euler-Lagrange equation 



(A2) 



so that for £^ < z < £^+\ the solution is completely deter- 



mined by the boundary conditions mi {t^, {£ a }) — 



and mi N) (i^+r, {£ a }) = 
write (for £^ < z < £^+1) 



Consequently we can 



( z; {£ a }) = F ^( z _ £^£^ +1 - (A3) 

for all fx and with some appropriate f unc tion which we 
need not determine. Substituting into (Al) then yields 



which is (4.19). The right hand side of (B4) is the inverse 
of the structure factor matrix element for a surface of 
fixed magnetization m x — (located at the centre of 
the thin film on the average) . If we calculate the binding 
potential W\ (£) for a single surface of fixed magnetization 
m x =0we can identify 



l(L) = W 1 (-) 



(B5) 



However from the stiffness-matrix formalism with N = 1 
we also find 



p0r = W{>(t) for£- 



dL 2 



(B6) 



valid in the soft mode phase T > T\y where z = 4. This 
completes the proof. 

A similar approach can be used for even systems and 
therefore we find generally 



K+ = 4 



d 2 l{L) 
dL 2 



(B7) 



JV+l 



W N {{£^})= Y^Wn-tn-l) 



exhibiting the required separable properties. 



(A4) 



APPENDIX B: CONNECTION WITH THE 
FREE-ENERGY 

Co nsider odd systems. We start by using ( |4.7D so that 
( 4.17 ) can be rewritten as 



rG (0,L) 



and using ( |4.7| ) again (with z = j) 
K. 



G (0,L) L 



G (04) 2 0V 2 



However from the 'product' relation (2.18) 



Go(0,|) 2 -Go(0,L)G o (^; 



so that 



K 4 



(Bl) 



(B2) 



(B3) 



(B4) 
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